Fluctuations and the Phase Transition in a Chiral Model with Polyakov Loops * 
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We explore the NJL model with Polyakov loops for a system of three colors and two flavors within 
the mean-field approximation, where both chiral symmetry and confinement are taken into account. 
We focus on the phase structure of the model and study the chiral and Polyakov loop susceptibilities. 
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I. INTRODUCTION 

Low-energy phenomena of QCD have been studied in 
various effective models based on chiral symmetry. How- 
ever, the relation of spontaneous chiral symmetry break- 
ing and confinement remains an open issue. Recently, 
color degrees of freedom were introduced in the Nambu- 
Jona-Lasinio (NJL) Lagrangian 0] through an effec- 
tive gluon potential expressed in terms of Polyakov loops 
(PNJL model) [i,l3- The basic ingredients of the model 
are constituent quarks and the Polyakov loop, which is 
an order parameter of the Z{3) symmetry of QCD in 
the heavy quark limit. The model has a non- vanishing 
coupling of the constituent quarks to the Polyakov loop 
and mimics confinement in the sense that only three- 
quark states contribute to the thermodynamics in the 
low-temperature phase. Hence, the PNJL model yields 
a better description of QCD thermodynamics near the 
phase transition than the NJL model. Furthermore, due 
to the symmetries of the Lagrangian, the model belongs 
to the same universality class as that expected for QCD. 
Thus, the model can be considered as a testing ground 
for the critical phenomena related to the breaking of the 
global Z{3) and chiral symmetries. 

It has been shown that the PNJL model, formulated 
at finite temperature and finite quark chemical poten- 
tial, well reproduces some of the thermodynamical ob- 
servables computed within lattice gauge theory (LGT). 
The properties of the equation of state [j| , the in-medium 
modification of meson masses Q as well as the validity 
and applicability of the Taylor expansion in quark chemi- 
cal potential used in LGT were recently addressed within 
the PNJL model SSBi]. In Ref. ^ the model was 
extended to a system with finite isospin chemical poten- 
tial and pion condensation was studied. 

Enhanced fluctuations are characteristic for phase 
transitions. Thus, the exploration of fluctuations is a 
promising tool for probing the phase structure of a sys- 
tem. The phase boundaries can be identified by the re- 
sponse of the fluctuations to changes in the thermody- 
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namic parameters. In this article we describe the suscep- 
tibilities of the order parameters and their properties in 
the PNJL model following Ref. 



II. NAMBU-JONA-LASINIO MODEL WITH 
POLYAKOV LOOPS 

An extension of the NJL Lagrangian by coupling the 
quarks to a uniform temporal background gauge field, 
which manifests itself entirely in the Polyakov loop, has 
been proposed to account for interactions with the color 
gauge field in effective chiral models 0, • The PNJL 
Lagrangian for three colors {Nc — 3) and two flavors 
(Nf = 2) with non-local four-fermion interactions is 
given by 
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where m = diag(m„, TO(j) is the current quark mass, /t = 
diag(/i„, ^d) is the quark chemical potential and t are the 
Pauli matrices. We assume isospin symmetry and take 
niu = md = mo and ^ t^^d ^ A*- 

Here we have introduced non-local interactions which 
are controlled by a form factor in order to deal with the 
ultraviolet singularities that appear in the loop integra- 
tions. In coordinate space, the form factor F{x) for the 
non-local current-current interaction reads: 
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A possible choice for the regulator is in momentum space 
given by [T^ : 
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where f{p) is the Fourier transform of the form factor 
F{x) and p is the three-momentum. The NJL sector 
is controlled by four parameters: the strength of four- 
fermion interaction Gs, the current quark mass mg and 
the constants a and A, which characterize the range of 
the non-locality. These parameters are determined in 
vacuum, for a given a, by requiring that the experimental 
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values of the pion parameters and the quark condensate 
are reproduced. 

The interaction between the effective gluon field and 
the quarks in the PNJL Lagrangian is implemented p.ip 
by means of a covariant derivative 
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where we introduce the standard notation = gA'^^ ^ . 
Here g is the color SU(3) gauge coupling constant and A° 
are the Gell-Mann matrices. 

The effective potential U of the gluon field in (j2.ip is 
expressed in terms of the traced Polyakov loop $ and its 
conjugate $ 



$ = — Tr^L . 
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where L is a matrix in color space related to the gauge 
field by 



L{x) = Vexp 
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with V being the path (Euclidean time) ordering, and 
(3 = 1/T with Ai = iAq. In the heavy quark mass 
limit, QCD has the Z(3) center symmetry which is spon- 
taneously broken in the high-temperature phase. The 
thermal expectation value of the Polyakov loop (<&) acts 
as an order parameter of the Z(3) symmetry. Conse- 
quently, ($) = at low temperatures in the confined 
phase and ($) ^ at high temperatures corresponding 
to the deconfined phase. 

The effective potential U{^, $) of the gluon field is 
expressed in terms of the Polyakov loops so as to preserve 



the Z{3) symmetry of the pure gauge theory [13 We 
adopt an effective potential of the following form 
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The coefficients Tq, and bi are fixed by requiring that 
the equation of state obtained in pure gauge theory on 
the lattice is reproduced. In particular, at To = 270 MeV 
the model reproduces the first order deconfinement phase 
transition of the pure gauge theory. 



III. SUSCEPTIBILITIES AND THE PHASE 
STRUCTURE 

In the PNJL model the constituent quarks and the 
Polyakov loops are effective fields related with the order 
parameters for the chiral and Z{3) symmetry breaking. 
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FIG. 1: The chiral Xmm (dashed-line) and the Polyakov loop 
Xii (solid-line) susceptibilities in the chiral limit as functions 
of temperature T for fj, = 0. 



In LGT the susceptibilities associated with these fields 
show clear signals of the phase transitions 14 1. 

In Fig. [1] we show the chiral Xmm and Polyakov loop 



Xii = {^^) — {^){^) susceptibilities computed at /i = 
in the PNJL model in the chiral limit within the mean 
field approximation. The chiral susceptibility exhibits a 
very narrow divergent peak at the chiral critical temper- 
ature Tch, while the peak of Xil is much broader and the 
susceptibility remains finite at all temperatures. This is 
due to the explicit breaking of the Z{3) symmetry by the 
presence of quark fields in the PNJL Lagrangian. Nev- 
ertheless, Xii still exhibits a peak structure that can be 
considered as a signal for the deconfinement transition in 
this model. One finds the interference of Xu with Xmm in 
Fig.[TJ At the chiral transition, T = Tch, there is a narrow 
structure in Xii- We stress that this feature is not related 
with the deconfinement transition, but expresses the cou- 
pling of quarks to the Polyakov loops. Thus, for the 
parameters used in the model, the deconfinement tran- 
sition, signaled by the broad bump in XiT, sets in earlier 
than the chiral transition at vanishing net quark density. 

The peak positions of the Xmm and XiT susceptibili- 
ties determine the phase boundaries in the (T, /^)-plane. 
At finite chemical potential, there is a shift of the chi- 
ral phase transition to lower temperatures. In Fig. [2] we 
show the resulting phase diagram for the PNJL model. 
The boundary lines of deconfinement and chiral symme- 
try restoration do not coincide. There is only one com- 
mon point in the phase diagram where the two transitions 
appear simultaneously. Recent LGT results both at van- 
ishing and at finite quark chemical potential show that 
deconfinement and chiral symmetry restoration appear 
in QCD along the same critical line [ij]. In general it is 
possible to choose the PNJL model parameters such that 
the both critical temperatures coincide at = 0. 

From Fig. [2] one finds that the slope of Tdec as a func- 
tion of ^ is almost flat, indicating that at low temperature 
the chiral phase transition should appear much earlier 
than deconfinement. However, there are general argu- 
ments, that the deconfinement transition should precede 
restoration of chiral symmetry (see e.g. p^. [l6j ) . In view 
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FIG. 2: The phase diagram of the PNJL model in the chiral 
limit. The solid (dashed) line denotes the chiral (deconfine- 
ment) phase transition respectively. The TCP (bold-point) is 
located at (T^ = 157, Hc = 266) MeV. 



of this, it seems unlikely that at T ~ the chiral sym- 
metry sets in at the lower baryon density than deconfine- 
ment. In the PNJL model, the parameters of the effec- 
tive gluon potential were fixed by fitting quenched LGT 
calculations. Consequently, the parameters are taken as 
independent on fi. However, it is conceivable that the 
effect of dynamical quarks can modify the coefficients 
of this potential, thus resulting in /x-dependence of the 
parameters. Consequently, the slope of Tdec as a func- 
tion of /i could be steeper ^. Consequently, the effective 
Polyakov loop potential (|2.7p should, with /^-independent 
coefficients, be considered as a good approximation only 
for ii/T < 1. 

While the susceptibilities Xmm and XiT exhibit ex- 
pected behaviors associated with the phase transitions, 
other Polyakov loop susceptibilities, 
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are negative in a broad temperature range above T^h 
(Tl| . This is in disagreement with recent lattice results, 
where \ii i s always positive in the presence of dynamical 
quarks [14| . A possible origin of this behavior could be 
the approximation to the effective Polyakov loop poten- 
tial used in the Eq. f^Tl} . 

Recently an improved effective potential with 
temperature-dependent coefficients has been sug- 



^ Such a modification was explored in Ref. [l7j where explicit fi- 
and Nf - dependence of Tq is extracted from the running coupUng 
constant as, using the argument based on the renormalization 
group. 
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FIG. 3: The xii = XiT susceptibility in the chiral limit as a 
function of temperature T for jj, — 0. The effective Polyakov 
loop potential (|3.2|l was used. 
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where 



a{T) ^ ao + ai 
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The polynomial in <I> and used in (|2.7|) . is replaced by 
a logarithmic term, which accounts for the Haar measure 
in the group integral. The parameters in (|3.2[) were fixed 
to reproduce the lattice results for pure gauge QCD ther- 
modynamics and for the behavior of the Polyakov loop. 
In Fig. [3] we show the xii susceptibihty calculated with 
the potential of Eq. (|3.2[) . The improved potential in- 
deed yields positive values for all the Polyakov loop sus- 
ceptibilities . We note that the phase diagram calculated 
with the improved potential is similar to that obtained 
with the previous choice of the Polyakov loop interac- 
tions, shown in Fig. [21 



IV. SUMMARY AND DISCUSSIONS 

We introduced susceptibilities related with the three 
different order parameters in the PNJL model, and ana- 
lyzed their properties and their behavior near the phase 
transitions. In particular, for the quark-antiquark and 
chiral density-density correlations we have discussed the 
interplay between the restoration of chiral symmetry and 
deconfinement. We observed that a coincidence of the 
deconfinement and chiral symmetry restoration is acci- 
dental. 

We found that, within the mean field approximation 
and with the polynomial form of an effective gluon poten- 
tial the correlations of the Polyakov loops in the quark- 
quark channel show an unphysical behavior, being neg- 
ative in a broad parameter range. This behavior was 
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traced back to the parameterization of the Polyakov loop 
potential. We argue that the Z(Af)-invariance of this po- 
tential and the fit to lattice thermodynamics in the pure 
gluon sector is not sufficient to provide correct descrip- 
tion of the Polyakov loop fiuctuations. Actually it was 
pointed out Q that the polynomial form used in this 
work does not possess the complete group structure of 
color SU(3) symmetry. The improved potential of Ref. Q 
yields a positive, i.e. physical, xu susceptibility, in qual- 
itative agreement with the LGT results. 
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